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ABSTRACT 

In order to determine if the material ablated from high- velocity clouds (HVCs) is a significant source 
of low-velocity high ions (C IV, N V, and O VI) such as those found in the Galactic halo, we simulate 
the hydrodynamics of the gas and the time-dependent ionization evolution of its carbon, nitrogen, and 
oxygen ions. Our suite of simulations examines the ablation of warm material from clouds of various 
sizes, densities, and velocities as they pass through the hot Galactic halo. The ablated material mixes 
with the environmental gas, producing an intermediate-temperature mixture that is rich in high ions 
and that slows to the speed of the surrounding gas. We find that the slow mixed material is a 
significant source of the low-velocity O VI that is observed in the halo, as it can account for at least 
~l/3 of the observed O VI column density. Hence, any complete model of the high ions in the halo 
should include the contribution to the O VI from ablated HVC material. However, such material is 
unlikely to be a major source of the observed C IV, presumably because the observed C IV is affected 
by photoionization, which our models do not include. We discuss a composite model that includes 
contributions from HVCs, supernova remnants, a cooling Galactic fountain, and photoionization by 
an external radiation field. By design, this model matches the observed O VI column density. This 
model can also account for most or all of the observed C IV, but only half of the observed N V. 

Keywords: Galaxy: halo — hydrodynamics — ISM: clouds — methods: numerical — ultraviolet: 
ISM 
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1. INTRODUCTION 

High ions from astrophysically abundant metals (e.g., 
C IV, N V, and O VI) in the interstellar medium 
(ISM) trace spatial and temporal transitions between the 
hot (T > 10 6 K) and warm/cool (T < few x 10 4 K) 
phases of the ISM. In the Galactic halo, above the 
disk, various physical structures may give rise to such 
transitio ns, including radiatiyely cooling Galactic foun- 
tain g as (lEdgar fc Chev alier 1986; Shap iro fc Benjamin! 
Il993t [Beniamin fc Shapirol I1993D . supernova remnants 
(|Sheltonl 119981 120061 ), evaporati ng cold clouds that are 
embedded in hot ambient g as (Bohring er fc HartquistJ 
Il987t iBorkowski et al.l Il990f ). or turbulent mixing lay- 
ers forme d where cool and ho t gas move relative t o one 
another (iBegelman fc Fabian! 119901: iSlavin et all 119931: 
lEsauivel et al.ll2006tlKwak fc Sheitonll2010D . Comparing 
the predictions of these various models with observations 
of high ions provides information on which physical pro- 
cesses are important in the Galaxy's halo. 

There are plenty of observational data for such com- 
parisons. High ions have been observed in the Galac- 
tic halo via their far-ultraviolet absorp tion lines, both 
at low velocities (5,90 km s" 1 ; e.g.. iPettini fc WestJ 
1982t ISavage fc Massal 119871: ISembach fc Savage! 1 1992 : 

[2 



Savage et all 119971. 120001 120031: iZsargo et al 



2003 



Indebetouw fc Shulll l200TlSavage fc Wakkeri 120091) and 
high velocities (>90 km s -1 : e.g.. | Murphy et alJ I 200 1 



Sembach etal]l2000L 
Collins et all 12004 



2003; 



2007 



?ox et all 12004 120051 1200C 
Shull et al.1 1201 lh . The 
low- velocity high ions' scale heights are »3-5 kpc 
(jSavage et all 19971 : iBowen et al.ll2008t ISavage fc Wakkeri 



l2009h . although the fi lling factor of the high-ion-rich 
material is small (e.g., IShelton et all 120071 ). C IV and 
O VI have also been observed in the halo via their emis- 



sion l ines (IShelton et al.ll200lLl2007ll2010HKorpela et al 



2006; Ottc & Dixon 2006; Dixon ct al. 2006; Welsh et al 
20071: iDixon fc SankritJl2008l : iPark et al.ll2009h .~ 

The observed ratios of high ion column densities arc of- 
ten used in attempts to di s tinguish betw e en models (e.g 



Sembach fc Savage] Il992t iSpitzerl 119961: ISembach et al l 
1997t ISavage et all 119971 : llndebetouw fc Shulll 12004) . 



leading to the conclusion that a single type of model 
cannot explain all of the observations. However, another 
important consideration is the quantities of each ion that 
the models predict. While the normalization of many of 
the aforementioned models is essentially a free parame- 
ter, it is important that the normalization required by 
the observations is physic ally reasonable (for example, 
llndebetouw fc Shulll [2004] point out that a worrisomely 
large number of turbulent mixing layers is needed to 
match the observed column densities). For other models, 
the normalization can be constrained in advance, with- 
out reference to the colu mn dens i ty me asurements. For 
example, in the model of lSheltonl (|2006f ). an ensemble of 
supernova remnants (SNRs) above \z\ — 130 pc can ac- 
count for 14-39% of the typical latitude-corrected O VI 
column densi t y obs erved toward extragalactic objects by 
ISavage eT al. (2003). In this case, the model's normaliza- 
tion is fixed using independently determined values for 
the rate and scale height of Galactic supernovae. 

Here, we consider a new model: high-velocity clouds 
(HVCs) interacting with their surroundings. HVCs 



2 



are interstellar clouds with |^lsr| > 90 km s 1 
([Wakker fc van Woerdenl fl99l . The first HVCs were 
discovered via their H I 21-cm emission (|Muller et al.l 
1 19631 ) , but ar e now know n also to have an ionized com- 
ponent (e.g., [Tuftl HOOl. On some sightlines, highly 
ionized high- velocity gas unassociated w i th hig h- velocity 
H I is observed (e.g., iSembach et~aH [2003J) HVCs 
may be material in a Galactic fountain (e.g., iBregmanl 
1980), material str ipped off satellite g alaxies (e.g., 
Gardiner fc Noguchil 11991 IPutman et aTT l2""004). mate- 
rial falling i nto the Galaxy fr om extragalactic space (e.g., 
IQort] 119661: iBlitz et al.l [19990 . o r material left over from 
the formation of the Galaxy ([Mailer fc Bullock! I2004Q . 
High ions can arise from the turbul ent mixing of cool 
cloud material with hot ambient gas (iSlavin et al.l 1 19931: 
Esauivel et~aTll2006t IKwak fc Sheltonl 120101: IKwak et all 



20m There are constraints on the number of HVCs in 
the Galactic halo. Given these constraints, we here con- 
sider whether or not material left behind by HVCs is a 
significant source of the low- velocity high ions observed 
in the halo. 

Our HVC simulations come from IKwak et al.l ([20TT[ 
hereafter Paper I), which used hydrodynamical sim- 
ulations to study the evolution of initially spherical 
HVCs traveling through a hot (10 6 K) ambient medium, 
thought to be representative of the upper halo. Our 
assuming the presence of this hot gas in the halo is 



back 
1995 


ground dBurrows & MendenhalU 119911: iWang & Yu 
: iPietz et al.l 119981: iWand 119981: iSnowdcn et al. 


1998 


, 2000; Kuntz & Snowdcn 


2000: Smith et al. 


2007 


: Galeazzi et al.l 2007; 


Hcnlcv & Shclton 


2008; 


Lei et al. 2009: Yoshino et al.l 


2009: iGuDta et al.l 


2009; 


Henlev et al. 2010) and of absorption lines, such as O VII 



and O VIII, in the X-ray spectra of active galactic nuclei 
(iNicastro et al. 2002; F ang et all 2003; Rasmusscn et all 



2003 1 IMcKernan et all 120041; I Williams et all 



Fang et al. 20061 : Bregman fc Llovd-Davie; 



2005 



2007 



Yao fc Wandl2007tlYao et alJl2008t 12009ft . The hot halo 



gas is likely due to Galactic fountains, with a possible 
contribution from accreting extragalactic material 
([Henlev et al.1 l20lol and references therein). However, 
it should be noted that the distribution of this gas is 
not currently well known, and it need not mostly fill the 
halo. 

In Paper I, we studied the high ions that arise when 
cool material ablates from an HVC and mixes with 
the hot ambient gas. To this end, we self-consistently 
traced the non-equilibrium ionization (NEI) evolution 
of carbon, nitrogen, and oxygen. In that paper, 
we concentrated on the high-velocity ions, those with 
line-of-sight speeds > 90 km s _1 . We showed that 
the column densities and column density ratios pre- 
dicted by our m odels overlapped those observed to- 
ward Comple x C (|Sembach et al.|[200l IPox et al.l 12001: 
iCollins et ai1l2007D . However, the high -ion-bearing ma- 
terial in our simulations eventually slows to ISM-like ve- 
locities, leading to copious quantities of low- velocity ions 
(with \v\ < 90 km s _1 ; see Section [2] for the exact defi- 
nition of "low velocity" used in this paper). In this pa- 
per, we describe the evolution of these ions. We use the 
Galactic HVC mass infall rate to estimate the column 
densities of high ions in the halo expected from HVCs in- 
teracting with hot ambient gas. When we compare these 



predictions with observations, we find that a significant 
fraction (>30%) of the observed low-velocity O VI could 
be due to this process. We therefore argue that any com- 
plete model of the high ions in the halo should take into 
account the contribution due to HVCs. 

The remainder of this paper is arranged as follows. In 
Section [5] we briefly describe our hydrodynamical model 
(see Paper I for more details). In Section [3] we describe 
the evolution of the low-velocity high ions, and how this 
is affected by our different model parameters. In Sec- 
tion [4] we present the average column densities of low- 
velocity high ions predicted by our models, and compare 
them with observations. In Section [5] we present O VI 
column density profiles predicted by our models. We 
discuss and summarize our results in Sections [6] and 
respectively. 

2. HYDRODYNAMICAL MODEL 

Our hydrody namical model is descr ibed in full in Pa- 
per I (see also IKwak fc Sheltonl [2010l) . Here, we give 
a brief overview of the model. The hydrodynamical 
simulations were carried out i n 2D cylindrical c oordi- 
nates using FLASH version 2.5 ([Frvxell et al.ll2000D . We 
tracked the ionization evolution of carbon, nitrogen, and 
oxygen using the FLASH NEI module^ The simula- 
tions include radiative cooling using the default FLASH 
cooling curve, which is a piecewis e pow er-law approxi- 
mation of the IRavmond fc Smith (119771) cooling curve 
(|Rosner et al.l 119781; IPeres et al.l ll982t T This cooling 
cur ye assumes collis i qnal i onization equilibrium (CIE); 
see IKwak fc Sheltonl (pOlOl) for some discussion of CIE 
versus NEI cooling rates. Note that the model does not 
include a magnetic field nor thermal conduction (see Sec- 
tion EXT]). 

In each simulation, the grid was initialized with a 
cool (10 3 K) HVC embedded in a hot (10 6 K) ambient 
medium. All of the clouds are initially spherical. The 
density of hydrogen atoms in each cloud (riH,ci in Ta- 
ble [lj pertains to the density of the cloud interior. Near 
the cloud's periphery, the density decreases smoothly 
to that of the ambient gas (nn.amb = ^H.ci/1000), and 
the temperature increases smoothly from 10 3 to 10 6 K, 
except in Model E, which has a sharp edge (see Fig- 
ure 1 in Paper I). Initially, the cloud interior, ambient 
medium, and transition zone are in pressure balance. 
Note that, although the densities are expressed in terms 
of the hydrogen number density, the gas also includes 
helium (uhc/^h = 0.1), so the total number density of 
atoms and ions is Linn- We assume that hydrogen and 
helium are fully ionized at all temperatures, so the elec- 
tron density n c — 1.2nn. 

The simulations were carried out in the initial rest- 
frame of the HVC; i.e., the HVC was initially at rest, 
while the ambient medium moved upward with velocity 
|v Zi ci|, where w ZiC i is the HVCs initial velocity in the ob- 
server's frame (the observer is assumed to be located be- 
low the domain, looking vertically upward). The bound- 
ary conditions allowed material to flow in from the bot- 
tom of the domain and flow off the top of the domain. 
Before extracting masses or column densities of high ions 

3 As noted in Paper I, we do not include the important high ion 
Si IV because it is more susceptible to photoionization. Modeling 
photoionization is beyond the scope of this paper. 
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Table 1 

Model Parameters 



Model 


ro 


fz,cl 








(pc) 


(km s _1 ) 


(cm -3 ) 


(Mq) 


(1) 


(2) 


(3) 


(4) 


(5) 


A 


20 


-100 


0.1 


120 


13 


150 


-100 


0.1 


4.9 x 10 4 


C 


150 


-150 


0.1 


4.9 x 10 4 


D 


150 


-300 


0.1 


4.9 x 10 4 


E 


150 


-150 


0.1 


4.9 x 10 4 


F 


300 


-100 


0.1 


4.0 x 10 5 


G 


150 


-100 


0.01 


4.9 x 10 3 


Note. 


— Column (1): 


Model identifiers. Col- 



umn (2): Initial cloud radius. For all except 
Model E, this radius is approximate, as the clouds 
do not have sharp edges (see Figure 1 in Paper I). 
Column (3): Initial velocity of the cloud along the 
2-direction measured in the observer's frame. Col- 
umn (4): Initial hydrogen number density of the 
cloud at its center. The ambient number density 
is 1/1000 this value. Column (5): Initial total 
mass of cloud, including hydrogen and helium. For 
all except Model E, this mass is approximate, as 
the clouds do not have sharp edges. Here, we use 
M[ n [ t T from Paper I — the initial mass of the ma- 
terial with T < 10 4 K. 

as a function of velocity, we transformed the velocities to 
the observer's frame. 

The parameters for each model in our suite of models 
are listed in Table Q] Model B is our reference model. 
The other models allow us to investigate the effects of 
cloud size (Models A, B, and F), cloud velocity (Mod- 
els B, C, and D), cloud density (Models B and G), and 
cloud density profile (Models C and E). 

In the following, we distinguish between high- and low- 
velocity ions using the same cuts in velocity that we used 
in Paper I. In Models A, B, F, and G (in which the initial 
velocity of the cloud was v ZyC \ = —100 km s -1 ), material 
with v z < —80 km s _1 is defined as high-velocity, while 
that with v z > —80 km s _1 is defined as low-velocity, 
where v z is the z-velocity of the gas in the observer's 
frame. For Models C, D, and E (in which w ZjC i was —150 
or —300 km s _1 ), the velocity cut is at v z = —100 km s -1 . 

When calculating the quantities of high ions that re- 
sult from th e clou d-ISM interactions, we assumed the 
IWilms et~a l. (2000[) interstellar abundances for carbon, 
nitrogen, and oxygen. Note that in Paper I we used 
older abundances, from I Allen] (|1973D . 

3. EVOLUTION OF LOW- VELOCITY HIGH IONS 

The hydrodynamical interaction between the model 
HVCs and the ambient gas is described in detail in Pa- 
per I. We begin this section by giving an overview of the 
processes. As a model HVC moves through the ISM, 
material is ablated from the cloud. The cool ablated 
material mixes with the hot ambient gas, creating gas of 
intermediate temperature (T ~ (1-3) x 10 5 K) in which 
high ions are abundant. The temperature of this mixed 
gas is affected both by radiative cooling and continued 
mixing with the hot ambient gas. The fractions of high 
ions in the mixed gas increase both by ionization of the 
initially cool ablated material, and by recombination of 
the initially hot ambient material. However, the fractions 
of all ions differ by varying degrees from those expected 
from CIE, as changes in the ionization balance lag be- 



hind the changes in the gas temperature that result from 
mixing and radiative cooling. Soon after material is torn 
from the clouds, it becomes relatively rich in high ions. 
At that time, the material is traveling almost as fast as 
the HVC, but as the ablated material continues to mix 
with the ambient medium, it slows. This causes the ve- 
locities of the ions to tend toward that of the ISM, and 
the ablated material to drift further from the cloud. 

In Figure [1] we show where the low- velocity ions are 
located relative to the HVC, by plotting the mass of 
each ion as a function of height. These plots were cre- 
ated using data from a modified version of Model B 
called Model B cxt , in which the domain extends up to 
•Zmax = 4400 pc, instead of z max = 2800 pc as in the orig- 
inal Model B (we increased z max for this purpose so we 
could trace the high ions over a greater range of heights). 
The height, z, is measured in the frame in which the HVC 
was initially at rest at z = 0. The HVC shifts upward 
during the course of the simulation, due to the ram pres- 
sure of the ambient medium. However, as this shift is 
only ~200 pc by t = 120 Myr, z is approximately the 
height above the HVCs current position. We show re- 
sults for two different epochs, t = 60 and 120 Myr. In 
each panel, the black curve shows the results of the NEI 
calculations, and gray curve shows the results obtained 
assuming CIE, where the ion fractions depend only on 
the local gas temperature. Note that, in general, there 
are significantly more high ions than are expected from 
CIE. Note also that the peaks seen at the two epochs do 
not represent the same ions - the peaks move through the 
domain at ~80-100 km s _1 , and so would move ~5-6 kpc 
in 60 Myr. 

We see that the low velocity C IV generally resides 
within 1400 pc of the cloud (for both CIE and NEI pre- 
dictions) at both 60 and 120 Myr. On average, both the 
N V and the O VI reside further behind the cloud. This 
is as expected - the mixed gas further behind the cloud 
contains a larger fraction of initially ambient gas and so, 
ignoring the effects of radiative cooling for the moment, 
is hotter and more highly ionized than that nearer to the 
cloud. 

Eventually, near the top of the domain, the mixed gas 
generally starts becoming too hot even for O VI, due to 
the continued mixing of the ablated cloud material with 
the hot ambient gas. This continued mixing and heating 
appears to be the ultimate fate for most of the ablated 
material - by the end of Model B, 38% of the material 
that was initially in the cloud and is now above z = 2 kpc 
is hotter than 6 x 10 5 K, while only 13% of this material 
is cooler than 1 x 10 5 K. The fraction of the initial cloud 
material that is hotter than 6 x 10 5 K by the end of 
Model B increases with increasing z. The hot mixed gas 
also becomes more tenuous, because of the low density 
of the ambient gas. Both of these effects lead to a fall off 
in the number of high ions. In addition, the mixed gas 
is closer to, although not yet in equilibrium. 

Although most of the mixed gas is hot, it was noted in 
Paper I that radiatively cooled mixed gas (T < 10 4 K) 
accumulates along the symmetry (r = 0) axis of the do- 
main. This radiatively cooled gas may be observable 
(e.g., via Ha emission), although such predictions are be- 
yond the scope of this paper. However, it should be noted 
that this accumulation of cooled gas along the symme- 
try axis may be an artifact of the cylindrical geometry. 



4 





1000 2000 3000 4000 

z (pc) 

Figure 1. Masses of low-velocity C IV, N V, and O VI (top to bottom) as functions of height, z, measured in the frame in which the 
HVC was initially at rest at z = 0. The data are from two epochs (60 and 120 Myr) of an extended version of Model B, Model B cx t, in 
which the top of the domain is at z = 4400 pc (instead of z = 2800 pc as in the original Model B). The mass of each ion was calculated by 
integrating over all cells at the specified height. The black and gray curves show NEI and CIE results, respectively. The data have been 
smoothed with a boxcar of full width 100 pc. 

Furthermore, because of the cylindrical geometry, this 
cooled gas represents a small fraction of the mixed gas 
mass. This also means that any further mixing of this 
cooled gas with hotter gas makes a negligible contribu- 
tion to the total high ion content. 

It is useful to calculate the time evolution of the num- 
bers of or masses of low-velocity high ions that result 
from the cloud-ISM interaction. Summing the amounts 
of such ions in the model domain at a given epoch pro- 
vides one measure, but the masses of the ions in the do- 
main are lower limits on the true masses of the ions that 
result from the HVC-ISM interaction. This is because 
these masses do not take into account material that has 
flowed out of the model domain (recall that each simu- 
lation is carried out in the cloud's rest frame and that 
material is allowed to flow off the domain at the top 
boundary) . 

For most models, we can make an estimate of the up- 
per limit on a given ion's mass by including material that 
has flowed off the domain. We are able to estimate the 
quantity of material that leaves the domain between each 
epoch in the model (i.e., the times at which the hydrody- 
namical data are output from the code) from the vertical 
velocity and position of the cells near the top of the do- 
main. Gas that is traveling upward at velocity v and 
that is a distance < vAt from the top of the domain at 
the current epoch will leave the grid by the next epoch, 
where At is the time between epochs (e.g., 0.5 Myr for 
Model B; note that At does not refer to the simulations! 



timestep, which is much shorter). We record the mass of 
each high ion that passes beyond the top of the domain 
at each epoch. 

Clearly, we cannot trace the evolution of the escap- 
ing material after it has left the computational domain, 
but we can tally the amount of a given high ion that 
has crossed the upper domain boundary at all preceding 
times. This tally serves as an upper limit on the amount 
of that ion beyond the upper boundary at any given time, 
because it ignores the possibility that some of the es- 
caped material further ionizes (as would occur when the 
gas mixes with additional hot ambient gas, raising the 
temperature above that which is favorable to C IV, N V, 
and O VI) or recombines. The sum of the amount of a 
given high ion that has moved out of the domain and the 
amount currently in the domain forms our upper limit 
on the total amount of that high ion present. 

Using Model B cxt , we investigated the evolution of 
the low-velocity high ions after they have risen above 
z = 2800 pc, the maximum height in Model B. We do 
this by following the time-evolution of the ion masses 
contained in select peaks in the Model B cxt mass-versus- 
height distribution (e.g., Figure [TJ as these peaks drift 
beyond the height of Model B. Figure [2] for example, 
shows the masses of O VI and C IV in a peak that 
we traced from t = 72 to 90 Myr, during which time 
it drifted from z w 2000 to 3600 pc. Both the mass of 
O VI and the mass of C IV in the peak decreased during 
this time (as did that of N V, not shown). Assuming that 
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Figure 2. Masses of low- velocity C IV (diamonds) and O VI 
(squares) in a single peak from the ion mass versus height distribu- 
tions (Figure [TJ as functions of time. The Model B GX t data (red) 
are from a peak that we traced from t = to = 72 to t = 90 Myr, 
during which time it moved from z 2000 to z ss 3600 pc. The 
ion masses include all ions within Az = ±250 pc of the peak. 
The Model C data (blue) are from a peak that we traced from 
t = to = 99 to t = 107 Myr, during which time it moved from 
2 S3 1300 to 2 fs 2500 pc. The ion masses include all ions within 
h.z = ±150 pc of the peak. 

this trend applies to all of the low- velocity high ions that 
rise past z — 2800 pc in Models B and B ext , the tally 
of all low-velocity high ions that have passed the upper 
boundary in Model B by some epoch of interest is an up- 
per limit on the mass of high ions above the boundary at 
that epoch. We expect this also to be true for the other 
models with v Z:C \ = —100 km s _1 (Models A, F, and G). 

In models with |v zc l| > 100 km s , the ram pressure 
of the ambient medium pushes the cloud upward during 
the course of the simulation (see Figure 4 in Paper I). 
As a result of this shift and the greater flow speed in the 
ambient medium, in the later stages of the simulation ab- 
lated material leaves the domain soon after it is torn from 
the HVC. The fractions of carbon, nitrogen, and oxygen 
in the C IV, N V, and O VI ionization levels, respectively, 
may still be increasing in this outflowing material when 
it leaves the domain. Also, in Model D in particular, a 
large fraction of the high ions that flow off the domain 
do so at high velocities. This high- velocity material may 
slow to low velocities after leaving the domain while re- 
maining rich in high ions. Thus, it is possible that the 
conclusion drawn from Model B cxt , above, docs not ap- 
ply to Models C, D, and E. We found, for example, by 
examining a Model C mass peak that the quantity of 
O VI ions was still increasing as the material approached 
the upper boundary of the domain (see Figure [5]), raising 
the possibility that it would have continued to increase 
with time after crossing the z — 2800 pc mark, had the 
simulation domain included greater heights. This was 
not the case for C IV or N V, whose masses had already 
begun to decrease while the mass peak was still in the 
domain (see Figure [2] for C IV). Hence, for Models C, 
D, and E, taking into account material that has left the 
domain may not always yield upper limits on the true 
total masses of the high ions. 

Figures [3] through [5] show the lower and upper limits 
on the masses of C IV, N V, and O VI, respectively, 
that result from the HVC-ISM interaction. These masses 
are plotted as functions of time for each of our seven 
hydrodynamical models. We plot the masses of both 
low-velocity and high-velocity ions (black and gray lines, 
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Figure 3. Mass of C IV as a function of time from our seven hy- 
drodynamical models (note the different ranges on the time axes). 
In each panel, the black lines show the masses of low-velocity ions, 
and the gray lines show the masses of high-velocity ions. The 
dashed lines show the mass of each ion that is in the model do- 
main, while the solid lines include ions that have escaped from the 
top of the domain (see text for details). Note that, in most panels, 
the gray solid and dashed lines are indistinguishable. 

respectively). In each case, we plot the mass of the ion in 
the domain (dashed line), and that mass plus the mass 
of the ion in the all material that has ever escaped from 
the domain (solid line). As noted previously, the former 
is the lower limit on the true mass, while the latter is our 
best estimate of the upper limit. 

3.1. General Features in the Evolution of the 
Low-velocity High Ions 

Here, we discuss the evolution of the low- velocity high 
ions in our reference model, Model B, and point out fea- 
tures that are common across all or most of the other 
models. In the following subsection, we will discuss how 
the different model parameters affect the evolution of 
these ions. 

The number of low-velocity high ions generally in- 
creases with time throughout the simulations. Admit- 
tedly, the masses of the high ions in the domain are de- 
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Figure 4. As Figure[3] but for N V. 

creasing toward the end of the Model B simulation. How- 
ever, as ^70% of the cloud's initial H I mass remains at 
the end of the simulation (Paper I), it is likely that this 
decrease is due to a relatively short-term fluctuation in 
the number of high ions, rather than due to a long-term 
decline in the processes that lead to high-ion-rich mate- 
rial. However, if we could run the simulation to much 
later times, we would eventually exhaust the cloud, and 
the processes that lead to high-ion-rich material would 
decline and eventually cease. The number of high ions 
would then decrease to zero as the high-ion-bearing gas 
either radiatively cooled or mixed with the million-degree 
ambient gas. 

When we compare the masses of low-velocity ions in 
Model B that do and do not take into account mate- 
rial that has escaped from the domain (solid and dashed 
black curves, respectively), we see that these lines di- 
verge by greater amounts as we go from C IV to N V to 
O VI (Figures El HJ and respectively). This behavior 
is understandable from Figure [TJ The low- velocity C IV 
mainly resides close to the HVC, well away from the top 
of the domain, and so very little low-velocity C IV flows 
off the domain. In contrast, the low-velocity O VI ex- 
tends further behind the HVC, nearer the top of the do- 



Model F 




150 
100 
50 



80 120 160 200 240 

time (Myr) 

Figure 5. As Figure[3] but for O VI. 

main, and so significant quantities of this ion escape from 
the domain. This general behavior is seen in the other 
models. However, it should be noted that in Model D 
the ambient material passes through the domain two or 
three times faster than in the other models. As a result, 
the mixed C IV-bearing gas is swept off the grid by the 
fast flow of the ambient gas. The fact that significant 
quantities of low- velocity C IV are lost from the domain 
in this simulation is reflected by the wide gap between 
the solid and dashed black curves in the Model D panel 
of Figure El 

We note that, in most cases, most of the high ions 
are at low velocities, despite their origin in the inter- 
action of an HVC with its surroundings. Although we 
are concentrating on the low- velocity ions here, we point 
out that the solid gray (high-velocity material includ- 
ing that which has left the domain) and dashed gray 
(high- velocity material in the domain) curves in Figures[3] 
throu.gb.0are indistinguishable in nearly every panel, im- 
plying that virtually no high-velocity high ions escape 
from the domain (i.e., the ions slowed to low velocities 
before they left the domain) Q Again, the exception to 

4 We remind the reader that "low-" and "high-velocity" refer 
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this is Model D. 

3.2. Effects of the Cloud Velocity, Profile, Density, and 

Size 

Cloud Velocity — Models B, C, and D have identical ini- 
tial clouds and ambient gas, but differing initial cloud 
velocities (—100, —150, and —300 km s _1 , respectively), 
allowing us to examine the effect of cloud velocity. As 
noted in Paper I, these velocities correspond to the sub- 
sonic, transonic, and supersonic regimes, respectively, 
resulting in significant differences in the morphological 
evolution of the HVCs. 

One key difference between the models is that in 
Model D a bow shock forms in front of the cloud, which 
helps protect the cloud from ablation at earlier times, 
and delays the onset of mixing (Paper I). As a result, 
there are very few high ions (high- or low-velocity) in 
Model D before t ~ 40 Myr. Another key difference is 
that the fast-moving ambient medium in Model D tends 
to sweep high ions off the domain before they are able to 
slow to low velocities in the observer's frame. As result, 
although high-velocity ions start becoming abundant in 
Model D at t ~ 40 Myr, low- velocity ions do not become 
abundant until t ~ 60 Myr. During this time delay, the 
HVC and hence the material ablated from it decelerated 
sufficiently for the mixed, high-ion-rich gas to reach low 
velocities before leaving the domain. 

After t ~ 60 Myr, despite the differences in the clouds' 
morphological evolution, the masses of low- velocity high 
ions that include the ions that have escaped from the 
domain (solid black curves in Figures [3] through [S]) agree 
within a factor of ^3 for all three models. If we con- 
sider only the ions in the domain (dashed black curves 
in Figures [3] through [S|) , we see that Model D has far 
fewer high ions than Model B after t ~ 80 Myr. This 
is likely because the Model D HVC shifts upward in the 
domain during the course of the simulation, due to the 
large ram pressure of the ambient medium (see Figure 4 
in Paper I). Eventually, the HVC gets so close to the top 
of the domain that, once again, the high ions are unable 
to slow to low velocities before leaving the domain. 

When considering the difference between the solid and 
dashed black curves in Figures [3] through [SJ Model C 
generally lies between Models B and D. As in Model D, 
the Model C cloud shifts upward in the domain during 
the simulation. Because the Model C cloud gets closer to 
the top of the domain, relatively more low- velocity high 
ions are swept off the top of domain than in Model B. 
This has the greatest effect on O VI, as O VI tends to 
exist further behind the HVC than C IV or N V. Model E, 
in which the cloud velocity is the same as in Model C, is 
similarly affected. 

Cloud Density Profile — The only difference between the 
initial parameters in Models C and E is the HVCs ini- 
tial density profile (smooth-edged versus sharp-edged; 
see Figure 1 in Paper I). The large density contrast at 
the edge of the Model E cloud inhibits the growth of 
shear instabilities and delays the onset of mixing, result- 
ing in fewer high ions than in Model C. However, after 

to velocities in the observer's frame, whereas the simulations were 
carried out in the HVCs initial rest frame. Hence, "low-velocity" 
ions are moving quickly in the simulation domain, whereas "high- 
velocity" ions move relatively slowly. 



t ~ 50 Myr, the predictions from these models that in- 
clude the ions that have escaped from the domain gen- 
erally agree within a factor of ~2, with Model E tending 
to yield larger masses. 

Cloud Density — The only differences between Models B 
and G are that the cloud and ambient densities arc 1/10 
as large in Model G than in Model B. As expected, the 
lower densities in Model G result in fewer high ions. How- 
ever, there are more high ions in Model G than one would 
expect from a simple rescaling of the Model B results. 
This difference is related to the fact that radiative cool- 
ing operates at a slower rate in Model G, because of the 
lower density. This difference in cooling rate affects the 
temperature distribution in the gas, which in turn affects 
the ionization and recombination rates. Nevertheless, de- 
spite these differences, the ion masses in Models B and 
G generally agree within a factor of ~2, after allowing 
for the factor of 10 difference in density. 

Cloud Size— Models A, B, and F differ by the HVCs 
initial radius (r = 20, 150, and 300 pc, respectively) and 
hence initial mass (Mhvc,o = 120, 4.9 x 10 4 , and 4.0 x 
10 5 Mq, respectively). Unsurprisingly, the more massive 
the HVC, the greater the mass of low- velocity ions that 
results from its interaction with its surroundings. 

We can examine this behavior more quantitatively by 
using a simple model of a uniform spherical HVC shed- 
ding mass at a rate proportional to its surface area (Pa- 
per I, Section 3.3.3). In this simple model, we expect the 
mass of a given low- velocity ion, M; on , at a rescaled sim- 
ulation time oit/ro, to be proportional to the cloud's ini- 
tial mass, Mhvc,o (see Equation (|A2[) in the Appendix). 
The rescaling of the simulation time by dividing by ro 
is necessary if we wish to compare models of different- 
sized HVCs at equivalent stages in their evolution, in 
the sense of their having had the same fraction of their 
initial H I mass ablated (see Paper I). Note that, in deriv- 
ing Mi OD (t/r Q ) oc Mhvc.Oj we have ignored the fact that 
the ions of interest will subsequently recombine or ionize. 
While this is a gross approximation, this simple model 
does provide some insight into the relative behaviors of 
Models A, B, and F. 

We find that the scaling of M- lon with Mhvc,o from our 
hydrodynamical models is close to the scaling expected 
from this simple model. When we compare Models B 
and A, we find M ion (t/r ) oc M^ vco , with 7 « 0.8-1.0 
(the exponent varies with t/ro and with ion). When we 
compare Models B and F, we find that 7 tends to be 
somewhat larger: 7 w 0.9-1.6. (In these comparisons, 
we used the ion masses that include material that has 
escaped from the domain, and we ignored early times in 
the simulations [t/ro < 0.2 Myr pc -1 , or t < 4, 30, and 
60 Myr for Models A, B, and F, respectively].) Thus, 
the values of Mi on (i/ro) scale approximately as expected 
with Mhvc.o- The deviations from Mi OD (t/r ) cx Mhvc,o 
may be due to our neglecting the recombination or ion- 
ization of the high ions. 

At later times in Model A, the masses of low-velocity 
N V and O VI that do and do not take into account 
material that has left the domain (black solid and black 
dashed lines, respectively) start to diverge significantly. 
This divergence is due to an effect already mentioned for 
Models C and E, above: the cloud shifts upward in the 



model domain, leading to more low- velocity high ions be- 
ing swept off the top of domain. In the case of Model A, it 
is the cloud's low mass-to-cross-section ratio that makes 
it susceptible to being pushed upward by the ambient 
medium's ram pressure. 

4. COMPARISON WITH OBSERVATIONS: THE HALO 
COLUMN DENSITY OF LOW- VELOCITY HIGH IONS 

4.1. Model Column Density Predictions 

In order to compare the results presented in the 
previous section with observations, we must calculate 
from these results the expected column densities of low- 
velocity high ions. In order to do this, we assume that the 
observed population of HVCs can be modeled by some 
number of model clouds, each of which is similar to one 
of our simulated clouds. In this case, 



A/hvc — 



■^HVC^HVC 
m HVC,0 



(1) 



is the number of model HVCs needed to account for the 
mass of high-velocity material, where Thvc and M^ivc o 
are the lifetime and initial H I mast0 of the model cloud, 

respectively, and .MhVc 1S tne observed infall rate of 
HVCs in the Galactic halo, expressed in terms of the 
total (neutral + ionized) hydrogen mass. By using an 
observed HVC infall rate that includes both neutral and 
ionized material, we are implicitly assuming that the ob- 
served HVCs began their lives as entirely neutral mate- 
rial, or, if they were initially partially ionized, that this 
does not have a significant effect on the quantities of high 
ions that result from their mixing with the hot halo. 

If (M ion ) H vc = J Mia n (t)dt/Tnvc is the time-averaged 
total mass of a given low- velocity ion due to ablation from 
a single model HVC, then the total mass of that ion in 
the halo due to ablation from the population of HVCs is 
given by 



M ion 

.halo — (M ion ) H vcM 



HVC 



•^HVC 
" i HVC,0 



M ion (t)dt. 



(2) 



If we assume that these ions are uniformly distributed in 
a cylindrical halo of radius Rmw above and below the 
disk, the average column density for a vertical sightline 
is 



iV(ion) = 



HVC 



2 7r -^MW^HVC,O r 



M ion (t)dt, (3) 



where raj on is the atomic mass of the ion (12.011, 
14.00674, and 15.9994 u for C IV, N V, and O VI, re- 
spectively [taken from the FLASH code]). 

The value of the HVC infall rate, A1h VC , is uncer- 
tain. If the observed HVCs are accreting extragalac- 
tic material, estimat es for -M^y^ ran ge from ~0.2 Mq 
yr- 1 (IMirabell H98l IPeek et al J 12008ft to ~1 M© yr" 1 
(|Wakker fc van Woerdenl Il997f ). If the observed HVCs 
are mainly due to a galactic fountain rather than due 
to accretion, there may be up to ~5 Mq yr" 1 circ ulat- 
ing through the halo dWakker fc van W ocrdcn 1997}). We 

5 Note that the cloud masses in Table[TJare total masses, and so 
must be divided by 1.4 to give the H I masses. 



assume A4^ vc =0.5 M© yr 1 . Considering that infall 
rates for individual co mplexes are as large as ~0.1 Mq 
yr- 1 for Complex C (jWakker et al.1 120071 : iThom et al.1 
EOQl) >0.1 Mq vr" 1 for t he S mith Cloud (usin g data 
from ILockman etall 120081 and IHill et al.1 I2009D . and 
~0.1-0.5 Mq yr -1 for the Magellanic Stream (the lower 
value is estimated from the mass of the negative - velocity 
portion of the St ream from iBriins et al .l 120051 : the up- 
per value is from IMirabell 11989ft . our assumed value of 

7W§ VC should b e reasonably c onservative. We also take 
R MW = 25 kpc (|Ferrierdl200l . 

Ideally, the time integration of M- lon (t) in Equation ^ 
should be carried out over the entire lifetime of the high 
ions in the halo. Note that this lifetime is, in general, not 
the time it would take for the HVC to be disrupted in the 
halo, because (a) the low-velocity high ions would tend 
to persist after the H I HVC has been completely dis- 
rupted, or (b) the HVC may reach the disk before it has 
been completely disrupted. In practice, if we integrate 
Mi on (t) up to the end of our simulations, then we obtain 
a lower limit on J Mi on (t)dt, as we are neglecting the 
contributions from beyond the ends of the simulations. 
In this subsection, we present the column density pre- 
dictions that result from the simulated part of the HVC 
evolution. In Section 14.31 below, we describe a method 
for including the estimated contributions from beyond 
the ends of the simulations. 

The column densities calculated using Equation ([3]) up 
to the times when the simulations ended are shown in the 
upper half of Tabled and are plotted in black in Figure^ 
In each case, we have calculated the column density using 
only the ions that are in the simulational domain ( "Do- 
main only" ) , and using the ions in the domain plus those 
that have escaped ( "Domain + Escaped" ) . For Model A, 
we also tabulate column densities resulting from only the 
first 16 Myr of the simulation (A* in Table [5J see Sec- 
tion l4.2[) . Note in particular that the Domain-only values 
are strict lower limits on the true model predictions (i.e., 
the column densities that the models would predict if 
we could trace all of the high ions indefinitely). This 
is because the mass of an ion within the simulation do- 
main is a lower limit on the true value of Mi on (i), and 
so each Domain-only column density is calculated from 
the lower limit of the integral of a quantity that is itself 
a lower limit. 

4.2. Model Comparison 

With few exceptions, the Domain-only and Domain + 
Escaped predictions for each model agree with each other 
within a factor of «2. The exceptions are all three ions in 
Model D, and O VI in Models A, C, and E. In these cases, 
the HVCs shift upwards, which allows large numbers of 
low-velocity high ions to escape from the domain during 
the course of the simulations. 

Comparing models having different velocities (Mod- 
els B, C, and D), we find that their Domain + Escaped 
predictions agree within a factor of 2, as do the Domain- 
only and the Domain + Escaped predictions for models 
having different density profiles (Models C and E). This 
is understandable, given the similarities in the ion masses 
between these models (see Section 13.2)) . The Domain- 
only and the Domain + Escaped predictions for models 
having different densities (Models B and G) also agree 
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Table 2 

Average Column Densities Predicted for the Galactic Halo 



N(C IV) JV(N V) N(0 VI) 

(10 13 cm" 2 ) (10 13 cm" 2 ) (10 13 cm' 2 ) 

Domain Domain + Domain Domain + Domain Domain + 

Model /9HVC a only Escaped only Escaped only Escaped 



A ■•■ 0.27 0.31 0.13 0.27 1.1 6.1 

A* ••• 0.088 0.088 0.043 0.044 0.35 0.41 

B ■•■ 0.26 0.30 0.10 0.13 0.99 1.5 

C ■•■ 0.14 0.18 0.043 0.072 0.32 0.96 

D ••• 0.031 0.17 0.015 0.076 0.11 0.91 

E •■■ 0.18 0.22 0.067 0.12 0.50 1.4 

F ■■■ 1.1 1.4 0.34 0.44 2.7 3.8 

G ■■■ 0.43 0.44 0.19 0.24 1.2 2.5 



Adjusted Column Density Predictions (see Section |4.3|I 



A 


0.721 


0.37 


0.43 


0.18 


0.37 


1.6 


8.5 


A* 


0.272* 


0.32 


0.32 


0.16 


0.16 


1.3 


1.5 


B 


0.216 


1.2 


1.4 


0.48 


0.60 


4.6 


6.8 


C 


0.344 


0.40 


0.51 


0.12 


0.21 


0.92 


2.8 


D 


0.257 


0.12 


0.65 


0.057 


0.29 


0.44 


3.6 


E 


0.228 


0.78 


0.97 


0.29 


0.54 


2.2 


6.1 


F 


0.327 


3.4 


4.3 


1.0 


1.3 


8,1 


12 


G 


0.272 


1.6 


1.6 


0.69 


0.87 


4.6 


9.3 



Note. — These column densities were calculated using Equation The column densities 
in the upper half of the table were calculated by integrating M; on (t) up to the end of each 
simulation, except where noted. As a result, these column densities are lower limits, because 
they neglect the contributions from ions at later times. The adjusted column densities in the 
lower half of the table are those from the upper half divided by /?hvc , to take i nto account the 
estimated contribution from ions beyond the ends of the simulations (see Section l4.3l for details). 
For each ion, the first column gives the predicted column density calculated only from the ions 
in the domain. The second column includes the ions that have escaped from the domain (see 
Section [3j . 

a The fraction of the HVC's initial H I mass that has ablated and/or ionized by the end of the 
simulation (except where noted), from Paper I. This quantity is not used in the upper half of 
the table, so we do not include it there. 

* Calculated by integrating Mi on (t) to t = 16 Myr, which corresponds to the same stage of 
evolution as the ends of Models B and F. Similarly, /3jjvc i s taken at t = 16 Myr, rather than 
at the end of the simulation. 



within a factor of 2 (note that, although M- lon (t) in Equa- 
tion ((3]) is lower for Model G, so too is AT|? y C ). 

Comparing the results from the different-sized model 
clouds (Models A, B, and F) is a little more compli- 
cated, as the simulated portions of the clou ds' lifetimes 
are not equivalent. As noted in Section 13.21 the simula- 
tion times must be divided by tq if we wish to compare 
different-sized HVCs at equivalent stages in their evolu- 
tion. Hence, we must reset the end time for the Mi on (i) 
integration, tf, such that tf/ro is the same for all clouds. 
We achieve this by using the full simulation times as tf 
for Models B and F (120 and 240 Myr, respectively), 
and setting tf to 16 Myr for Model A. The column den- 
sity predictions resulting from the new integration of the 
Model A results are labeled A* in Table [2J 

Having obtained column density predictions for equiv- 
alent portions of the clouds' lifetimes, we can compare 
them to the expectations for the simple model of a spher- 
ical cloud mentioned in Section[3T2] In this simple model, 
if tf/rQ is the same for all clouds, we expect iV(ion) oc ro 
(see Equation (|A3I) in the Appendix). The ro ratios for 
Models A*:B:F are 0.13:1:2 (Table [T]). In comparison, 
the ratios of the Domain + Escaped column densities 
are 0.29:1:4.7 for C IV, 0.34:1:3.4 for N V, and 0.28:1:2.6 
for O VI. The ratios of the Domain-only column den- 
sities are generally similar: 0.34:1:4.3, 0.41:1:3.3, and 
0.35:1:2.8 for C IV, N V, and O VI, respectively. Hence, 
the column densities increase with rg, but not linearly. 



The fact that the column densities for Models A* and F 
are somewhat larger relative to the Model B values than 
expected from tq is connected to the observation that the 
relationship between Mj on (t/ro) and M^y C is generally 
somewhat shallower than linear proportionality between 
Models A* and B (7 « 0.8-1.0), and generally some- 
what steeper than linear proportionality between Mod- 
els B and F (7 ss 0.9-1.6; see Section EHj) . 

For Model A, we can extend the integral of M ion (t) to 
later times, and naturally iV(ion) increases as a result. 
However, we do not have data from Models B and F for 
later times with which to compare the Model A results. 

4.3. Including the Contribution from Ions Beyond the 
Ends of the Simulations 

Here, we describe a method for estimating the con- 
tributions to the column densities from ions after the 
simulations have ended. From Equation ([3]), we see 
that the predicted column densities depend on the in- 
tegral J M ion (t)dt. Ideally, the time integration should 
be carried out over the entire lifetime of the high ions 
in the halo, but in Section 14.11 we carried out the in- 
tegration only up to the end of each simulation at time 
t = T slm . To estimate the contribution from ions be- 
yond the ends of the simulations, we assume that each 
and every unit mass of H I "lost" from the HVC to 
ablation and ionization makes an approximately equal 
contribution to J Mi on (t)dt. In this case, the integral 
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Figure 6. Comparison of the observed and predicted low-velocity 
halo column densities for C IV, N V, and O VI (top to bottom). 
The lower and upper ends of the black vertical bars show the "Do- 
main only" and "Domain + Escaped" predictions, respectively, for 
each model from the upper half of Table [2] For Model A, we plot 
the value obtained by integrating M- lon (t) up to the end of the sim- 
ulation (i.e., A from Table [2] rather than A*). These predicted 
column densities are lower limits, as we did not follow the high 
ions for the entire lifetime of each cloud. The gray vertical bars, 
offset slightly to the right, show the adjusted column density pre- 
dictions from the lower half of Table [2] These predictions take into 
account the estimated contr ibut ion from ions beyond the ends of 
the simulations (see Section I4.3I for details). The horizontal solid 
lines show the observed halo column densities. For C IV we plot 
the average value of log[iV(C IV) sin|b|1, with the dashed lines in- 
dicating the standard deviation (Savage & Wakker 2009, Table 4, 
specifically the values for all extragalactic sightlines). For N V 
we plot the best-fit value of the midplane ion density multiplied 
by the io n scale height, with the dashed lines indicating the un- 
certainty (Savage ct al. 1997, Table 7, specifically the values that 
properly takes into account the fact that some sightlines yielded 
only upper limits on the N V column density). We assume that 
the contribution from the Local Bubble is negligible for these ions. 
For O VI we plot the average value of log[7V(0 VI) sin |6|], with the 
dashed lines indicating the standard deviation ( Savage ct al~1 l2003 . 
Table 3, specifically the values for the full dataset). We have sub- 
tracted 7 X 10 12 cm~ 2 to all ow for the contribution from the Local 
Bubble (Oegerle et al. 20Ql). 



be replaced by J Tsim M\ 



,lost 



A^HVC,o) 



J Tsim Mi on (t)dt/0HVc, where M^y C lost is the mass of H I 
lost from the HVC to ablation and ionization, and /?hvc 
is the ratio of the lost mass to the HVC's initial H I mass 
(defined in Paper I). Substituting this new integral into 
Equation ([3]) allows us to approximate the column den- 
sities that would be predicted if we could follow the full 
evolution of the cloud. In practice, this substitution can 
be accomplished by dividing the column densities calcu- 
lated in Section |4~T1 by /?hvc- 

The lower half of Table [5] contains the adjusted column 
density predictions for Models A through G, obtained by 
dividing the values in the upper half of the table by /3hvc ■ 
These adjusted column density predictions are also plot- 
ted in Figure [6l in gray. The adjustment increases the 
column densities for Models B through G by factors of 
~3-5 (-0.5-0.7 dex), and those for Model A by «40% 
(0.14 dex). The lower half of Table [5] also contains the 
adjusted column densities for Model A*, obtained by di- 
viding the relevant values in the upper half of the table 
by the value of /3hvc at i = 16 Myr. 

After dividing by /?hvCj the Domain-only column den- 
sity predictions for Model A are systematically higher 
than those for Model A*. As the Model A predictions 
were obtained by integrating the ion masses to a later 
time than the Model A* predictions, this discrepancy 
suggests that we are still underestimating the contribu- 
tion from ions at times later than the ends of the sim- 
ulations. Nevertheless, the Domain-only predictions for 
Models A and A* agree within 25%, which gives us some 
confidence that dividing by /3hvc leads to reasonable es- 
timates of the contribution of ions from beyond the ends 
of the simulations to the column densities. It should be 
noted that the Domain + Escaped predictions for Mod- 
els A and A* in the lower half of Table [5] are less con- 
sistent than the Domain-only predictions, especially for 
N V and O VI. However, the Domain + Escaped predic- 
tions for Model A may be significantly overestimated, as 
relatively large quantities of low-velocity N V and O VI 
were swept off the top of the domain in this simulation, 
and we were unable to follow these ions' subsequent evo- 
lution. Hence, the discrepancy in the adjusted Domain 
+ Escaped predictions for Models A and A* does not 
argue against the reliability of our method for estimat- 
ing the contribution from ions beyond the ends of the 
simulations. 

It should be noted that there is a time delay between a 
unit mass of cool material being stripped from the cloud 
and that material becoming rich in high ions. As a re- 
sult, the above-described division by /?hvc will neglect 
the contribution to J M lon (t)dt from material that has 
been stripped from the cloud, but has not yet completed 
its ionization evolution by the end of the simulation. 
However, as this neglected contribution to J Mi on (t)dt 
is likely to be small, we do not attempt to further adjust 
our column density predictions to take it into account. 

4.4. Comparison with Observations 

Owing to the large number of models examined and 
the multiple methods for estimating their contributions 
to the column densities of low-velocity high ions, our 
column density predictions cover a wide range. At the 
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high end, they can account for > 1/3 of the observed O VI 
column density, while our lower limits account for only a 
few percent of the observed C IV. 

Let us first consider the column densities from the up- 
per half of Table [5] (the values that exclude the contribu- 
tions from the times beyond the ends of the simulations) . 
These are plotted in black in Figure [5J where they are 
compared with the observed values (10 x 10 13 , 2.5 x 10 13 , 
and 1 5.5 x 10 13 cm" 2 for C IV, N V, and O VI, respec - 
tively [Savage fc Wakker1[200l ISavaee et al.lll997l [2003j : 
see figure caption for details). In general, we find that 
these model column densities account for only a small 
fraction of the low-velocity C IV, N V, and O VI ob- 
served in the halo. The exceptions are Models A and F, 
whose Domain + Escaped O VI predictions agrees with 
the observed value within a factor of ^2.5 and ~4, re- 
spectively. 

When we include the estimated contribution from the 
times beyond the ends of the simulations (tabulated in 
lower half of Tabled plotted in gray in Figure [5]), we see 
that several models' predictions come within a factor of 
~3 of the observed O VI column density. The Domain- 
only and Domain + Escaped predictions from Models B 
and G are within a factor of 3.5 of the observed value, 
while the Domain-only and Domain + Escaped predic- 
tions from Model F and the Domain + Escaped predic- 
tion from Model A are within a factor of 2 of the ob- 
served value. The Domain-only predictions are of partic- 
ular interest, as they are lower limits on the true model 
predictions. This is because the original Domain-only 
predictions are strict lower limits on the true predictions 
(see Section |47T|) . and dividing by /?hvc may underesti- 
mate the contribution from times beyond the ends of the 
simulations (see Section |4~3")) . Hence, our HVC models 
can account for a significant fraction of the low-velocity 
O VI observed in the halo (e.g., our reference model, 
Model B, can account for 30-44% of the observed O VI, 
where the lower limit is taken from the Domain-only pre- 
diction, and the upper limit is taken from the Domain + 
Escaped prediction). 

For C IV and N V, we find that fewer of our models 
predict column densities that are within a factor of ^3 
of the observed values. The Domain + Escaped C IV 
and N V predictions from Model F that include the con- 
tributions from times beyond the end of the simulation 
agree with the observed values within factors of 3 and 2, 
respectively, while the equivalent Model G N V predic- 
tion agrees with the observed value within a factor of 3. 
The other models typically account for <15% and <25% 
of the observed C IV and N V, respectively (e.g., our 
reference model, Model B, can account for only 12-14% 
of the observed C IV, where the lower and upper limits 
are again taken from the Domain-only and Domain + 
Escaped predictions, respectively). We will discuss these 
results further in Section |6] 

5. COLUMN DENSITY PROFILES 

In this section we present O VI column density pro- 
files (i.e., column density as a function of velocity) de- 
rived from one of our HVC models, and compare them in 
general terms with the observed profiles. Note that, al- 
though the main topic of this paper is low-velocity high 
ions, these column density profiles span both low and 
high velocities. 



Figure [JJ shows O VI profiles for four different ver- 
tical sightlines from two epochs of Model B, and Fig- 
ure [8] shows profiles for five different horizontal sightlines 
from one epoch of the same model. The column den- 
sity of each hydrodynamical cell along the line of sight 
was calculated from the gas density and the relevant ion 
fraction in that cell; the velocity bin was chosen accord- 
ing to the cell's line-of-sight velocity, v, in the observer's 
frame (for the vertical and horizontal sightlines, the ob- 
server is located below and to the right of the model 
domain, respectively). The dotted column density pro- 
files in Figures [JJ and [5] do not include thermal broaden- 
ing, and were constructed by summing the contributions 
from all the cells along the line of sight, assuming that 
each cell's column density profile is a Dirac 8 function at 
v. The solid column density profiles do include thermal 
broadening, which we approximated by convolving each 
cell's column density profile with a boxcar of full width 
26, where b = (0.129 km &- l )y/(T/K)/A is the thermal 
velocity-spread parameter, T is the gas temperature in 
the cell, and A is the r elative atomic mass of the ion in 
question (|SpitzeH ll978. Equation 3-21). Note that the 
profiles in Figures [7] and [5] are much bro ader than the 
FUSE spectral resoluti on (~20 km s _1 ; iSahnow et al.1 
[2lMlWood"eral1[2l)0l . 

Figures [7J and [5] show that the column density profiles 
with and without thermal broadening are similar to each 
other, regardless of epoch or sightline. This means that 
the broadening of the profiles is due to variation in the 
bulk fluid velocity along the line of sight, rather than 
thermal broadening. For the horizontal sightlines, this 
variation in the bulk fluid velocity is due to turbulence 
in the mixed gas. For the vertical sightlines, however, 
there is an additional effect. The ions' velocities tend to- 
ward that of the ambient medium (v = in the observer's 
frame) the further they are from the HVC. Therefore, a 
vertical sightline through our model domains will typi- 
cally sample both high- and low-velocity ions, spanning 
a continuous range of velocities. As a result of this addi- 
tional effect, the column density profiles for the vertical 
sightlines are generally somewhat broader than those for 
the horizontal sightlines. 

FUSE observations of O VI absorption in the halo 
show some qualitative agreement and some qualitative 
disagreement with our predicted column density profiles. 
The observed line widths are generally much broader 
than the expected thermal width, as in our model pro- 
files. The observed velocity-spread parameters for low- 
and high-velocity O VI absorption ar e typically b ~ 50- 
70 and ~30-50 km s -1 , respectively (|Savage et al.H2003l 
Scmbach ct al. 2003), compared with b = 18 km s _1 for 
thermal broadening at T — 3 x 10 5 K. For several sight- 
lines, absorption is observed over a continuous range of 
velocities from y ~ to v ~ — 200 or +200 km s _1 
(see Figure 1 in lWakker et aLlfeOOSt ). in qualitative agree- 
ment with our model column density profiles for vertical 
sightlines. However, the observed column density profiles 
typically have relatively more ions near v = than our 
model profiles. These extra ions near v — are likely due 
to other sources of low-velocity high ions, in addition to 
those predicted by our HVC models. We discuss further 
the idea of multiple sources of high ions in the halo in 
Section O 
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Figure 7. O VI column density profiles taken from two epochs (left: 60 Myr, right: 120 Myr) of Model B. The column density profiles 
are for material lying along vertical sightlines at four different impact parameters relative to the cloud center (top to bottom: 10, 80, 160, 
and 320 pc). The solid and dotted lines show the profiles with and without thermal broadening (see text for details). 



6. DISCUSSION 

We have analyzed existing simulations of HVCs (Pa- 
per I) in order to estimate the quantities of low-velocity 
ions that result from the passage of HVCs through the 
hot halo. In Section 14.41 we showed that some of our 
models could account for a significant fraction of the 
observed low-velocity O VI, but not of the low-velocity 
C IV or N V. In this section, we first discuss our model 
assumptions, in particular neglected physical processes 
(Section |6.1.1[) , the pressure and temperature of the halo 
(Section l6.1.2p , and the factors in Equation ([3]) for which 
we had to assume values (Section 16.1.31) . Then, in Sec- 
tion H21 we discuss our model alongside other models of 
the high ions in the halo. 

6.1. Model Assumptions 
6.1.1. Neglected Physical Processes 

Our model does not include a magnetic field, as our 
2D geometry prevents our modeling realistic field con- 
figurations. Magnetic fields are known to suppr ess the 
development of turbulence (e.g.. lRvu et alJfeOOOf ). while 
the development of turb ulence differs in 2D and 3D (e.g., 
iStone fc Normanl lT992l). It is unclear whether adding a 
magnetic field and a third dimension to the model would 
result in more or fewer high ions, compa red to our cur- 
rent simulations. iKwak fc Sheltonl (|2010t ) point out that 
their 2D hydrodynamic simulations of plane-parallel mix- 
ing layers predict ion column density ratios that are sim- 
ilar to those pr edicted by the 3 D ma gnetohydrodynamic 
simulations of Es quivel et al.l ((2006) . However, we are 
unable to make a similar statement regarding the mag- 
nitudes of the column densities. 



Our model also does not include thermal conduction. 
Like a magnetic field, thermal conduction may sup - 
press the growth of instabilities ([Orlando et al.1 12008). 
although we cannot predict the extent to which thermal 
conduction would affect the development of turbulence 
in our simulations. Assuming that turbulence is able to 
develop, the inclusion of thermal conduction in our simu- 
lations would be unlikely to greatly affect our results, as 
turbulent diffusion of heat should dominate over diffusion 
by thermal conduction (see Section 5 of Paper I). 

6.1.2. The Pressure and Temperature of the Halo 

Our model assumes that the halo pressure is uniform 
over the distance that the cloud falls through the halo 
(>10 kpc, assuming a cloud speed of 100 km s _1 and 
a lifetime of >100 Myr). This would be an unrealistic 
assumption for HVCs traveling through the lower halo, 
where the total interstellar pressure more than doubles 
from 3600 to 8 700 cnT 3 K between \z\ = 3 and 1 kpc 
(|Ferrierd lT998). However, the ambient halo pressure in 
our model is much lower than this (230 cm" 3 K in most 
models, 23 cm -3 K in Model G), and may be more ap- 
propriate for HVCs traveling through the upper halo 
(\z\ > 10 kpc). Although the pressure in the upper halo 
is uncertain, in Paper I we noted that our chosen ambi- 
ent density (n H = 1.0 x 10~ 4 cm~ 3 for Models A-F) is 
similar to some previous observational estimates of the 
density in the up per halo (~(l-few) x 10~ 4 cm" 3 for 
\z\ > 10 kpc: e.g.. IWemer fc WillianTsI [19961 iPeek et al.l 
I2007t iGrcevich fc Putmanll2009D . As these observational 
estimates of the density in the upper halo vary by only a 
factor of a few over a wide range of heights above the disk 
(tens of kpc) , the pressure gradients will not be large in 
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Figure 8. As Figure [7] but for horizontal sightlines through the 
Model B domain at t = 120 Myr. Profiles are plotted for sightlines 
at five different heights, measured relative to the cloud's initial po- 
sition in the domain (top to bottom: 0, 240, 410, 730, and 1800 pc). 



the upper halo if the temperature is reasonably uniform 
in this region. 

The orig i n of HVCs remains uncertain, although 
I Peek et alT (|2007|) point out that the observation of 
HVCs >10 kpc above the disk argues aga inst their ori- 
gin in a Galactic fountain (Bregman! 1980[ ). Some HVCs 
have low metallici ties (e.g., ~0.13 solar for Complex C; 
iCollins et~a l. 2007), which also argues against their being 
composed primarily of fountain material. If the major- 
ity of HVCs originate in extragalactic or circumgalac- 
tic material, they would travel great distances through 
the upper halo, in a relatively low pressure, low density 
environment similar to that in our models. As noted 
above, the pressure gradients may not be large in the 
upper halo. Furthermore, clouds that have transverse 
components to their motion through the halo would ex- 
perience less change in ambient pressure in a vertically 
stratified halo than if they were traveling vertically down- 
ward. Modeling the change in ambient pressure as an 
HVC falls into the Galaxy is beyond our current models, 
but could be incorporated into future simulations. 

Our model also assumes a halo temperature of 1 x 
10 6 K. X-ray observations indicate the presence of 
^million-degree gas in the halo, although its tempera- 
ture structure and filling factor remain uncertain. We 
did not investigate different ambient temperatures in 



our suite of models. However, in their simulat i ons o f 
2D plane-parallel mixing layers. iKwak fc Shelter] (f2010h 
found that increasing the temperature on the hot side 
of the interface from 1 x 10 6 to 3 x 10 6 K did not have 
a large effect on the ion column densities for sightlines 
p erpendicular to th e inter face. 

IKwak k. Shelter] (|2010fl did not carry out a corre- 
sponding simulation with a temperature lower than 1 x 
10 6 K. However, if the halo temperature were much lower 
than 1 x 10 6 K, the halo's ROSAT R2/R1 ratio would 
be lower than the observed value (e.g., R2/R1 = 0.48 
for a 5 x 10 5 K plasma, while observational analyses 
yield R2/R1 > 0.5 for the halo's intrinsic emission; 
ISnowden et aTlEfflOOl IKuntz fc Snowdei] [20001) . If the 
ambient temperature were even lower, then the ambi- 
ent medium would not persist in the upper halo, as 
it would have a relatively short cooling time. For ex- 
ample, for an ambient medium with T — 3 x 10 5 K 
and riH = 1-0 x 10~ 4 cm~ 3 , the cooling time would 
be ^100 Myr (calculated using the 1993 version of the 
iRaymond fc Smithl 119771 code). This time is similar to 
the lengths of our simulations. 

6.1.3. -Mhvc> Rmw, and Elemental Abundances 

When we used Equation ([3]) to calculate the column 
density of a given ion from our hydrodynamical simula- 
tions, we had to assume values for three important fac- 
tors: Al^vc, Rmw, and the elemental abundance (the 
abundance affects the values of M- lon {t) derived from our 

simulations). As noted in Section l4~T| the value of -M§ vc 
is uncertain. However, our assumed value of 0.5 Mq 
yr _1 should be reasonably conservative (in the sense of 
tending to lead to underestimates of the predicted col- 
umn densities), given the large infall rates associated 
with individual complexes (e.g., Complex C, the Smith 
Cloud, and the Magellanic Stream). The radius Rmw 
over which the halo ions arc distributed is more uncer- 
tain. By choosing Rmw = 25 kpc, we have assumed 
the ions that result from HVCs interacting with the hot 
halo are spread uniformly above the entire stellar disk. 
However, if these ions tend to be concentrated toward 
the center of the Galaxy, which is a reasonable assump- 
tion, then our choice of Rmw should also be reasonably 
conservative. 

For the e lemental abundance s, we used the interstellar 
values from lWilms et al.l ()2000l ) , which are in good agree- 
ment with r ecent measurements of solar photo spheric 
abundances IlLodderdfeOOa lAsplund et al.l 120091 ) Fl The 
iWilms et alJ (|2000| ) interstellar abundances of carbon 
and oxygen include atoms that are in dust: the gas-phase 
abundances of carbon and oxygen in the local ISM are 
0.23 and 0.1 9 dex lower than the IWilms et al.l va lues, 
respectively (|Cardelli et"aT1ll996tlMever et al.lll998h. Ni- 
troge n, however, is not depleted onto dust ( Mever et alJ 
1997). Carbon, nitrogen, and oxygen abundances in the 
halo are more uncertain, but in the lower halo at least 
(|z| < 2 kpc), total abundances appear to be approxi- 
mately solar (based on the gas-phase sulfur abundance), 
and there is less depletion onto dust than in the disk for 

6 Note that th e ion masses and column densities in Paper I were 
calculated using Allen ( 1973) abundances , which are 0.14, 0.08, 
and 0.13 dex larger than the Wilms et al. (2000) abundances for 
carbon, nitrogen, and oxygen, respectively. 
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elements heavier than oxygen ([Savage fc Sembachlll996l ). 
If the carbon, nitrogen, and oxygen abundances are also 
approximately solar in the halo, and less depleted onto 
dust than in the disk, then our M- lon {t) predictions for 
C IV and O VI would need to be revised downward by 
no more than ~40% to compensate for the possibility 
that these elements are depleted onto dust, while our 
N V predictions would not need revision. 

However, it should be noted that varying the abun- 
dances would affect the cooling curve used in our sim- 
ulations. While the FLASH manual is not explicit, the 
FLASH coolin g curve is lik ely based on a cooling curve 
calculated with I Allen! (1 19731) abunda n ces, a s these are the 
defaults for the iRavmond fc Smith! (|1977| ) code. Mod- 
ifying the cooling curve would affect the quantities of 
high ions predicted by our simulations, by affecting the 
amount of gas at optimal temperatures for such ions. 
Quantifying the extent of this effect is beyond the scope 
of this paper. However, while lowering the abundances 
would reduce Mi on (t), it would also reduce the cooling 
rate, meaning that gas would remain at high-ion-rich 
temperatures for longer. As a result, lowering the abun- 
dances would not necessarily lead to a commensurate 
reduction in J AL ion (t)dt, and hence in the predicted col- 
umn densities. 



A Composite Model of the High Ions in the 
Galactic Halo 



z\ , which appears not to be the case observa- 



We pointed out in Section 14.41 that the Domain-only 
column density predictions in the lower half of Table [2] 
are lower limits on the true model predictions. In the 
previous subsection, we argued that our choices of -M^yc 
and Rmw should be reasonably conservative, and that 
adjusting the elemental abundances is unlikely to result 
in large adjustments in the predicted column densities. 
Hence, our statement in Section POl that our HVC models 
can account for a significant fraction of the low-velocity 
O VI observed in the halo is quite robust. 

While the Domain-only predictions from the lower half 
of Table [2] are lower limits, the corresponding Domain + 
Escaped predictions are generally upper limits on the ex- 
pected ion column densities. This is because, for most 
models, the numbers of low-velocity ions in the mate- 
rial that has flowed off the domain decreases with time 
(see Section [3]). Therefore, taking these lower and upper 
limits for Model B (our reference model), we find that 
HVCs could account for 30-44% of the low-velocity O VI 
in the halo, but only 12-14% of the low-velocity C IV 
(Section 

The preceding statement raises the question, why can 
our model account for a significant fraction of the O VI 
in the halo, but not of the C IV? The answer is likely be- 
cause our model does not include photoionization, which 
can increase the amount of C IV relative to O VI, due 
to C IV's lower ionization potenti al. For example, the 
coolin g Galactic fountain model of iShapiro fe Benjamin! 
(1993), which includes photoionization, predicts ~2-7 
times as muc h C IV relative to O VI as the cooling foun- 
tain model of lEdgar fc Chevalierl (|1986l ) , w hich does not. 
Mode ls that include photoionization (e.g.. Hto fc Ikeuchil 
1988) can explain the enhancement in C IV (and Si IV, 
which is not studied in this paper) rela tive to N V ob- 
served at large \z\ (jSavage et al.l Il997h . However, it 



should be noted that the llto fc Ikeuchil (|1988) model also 
predicts that C IV should be enhanced relative to O VI 
at large \z\, which appears i 
tionally (Sa vage et al.M2003f ). 

As our HVC model does not include photoionization, 
and therefore tends to significantly underpredict the 
amount of C IV in the halo, we are not putting it for- 
ward as the only source of high ions in the halo and as 
an alternative to other models. Instead, we are pointing 
out that a complete model of the Galactic halo should in- 
clude the contribution of HVCs to the low-velocity high 
ions, particularly O VI. While developing a complete self- 
consistent model of the high ions in the halo is beyond 
the scope of this paper, we will give an example of how 
different existing models of the halo high ions could be 
pieced together. Note that we consider only the total 
column densities of the ions, not their z distributions, as 
we cannot derive the z distributions of the ions from our 
current HVC model. However, if the HVCs spend most of 
their lifetimes in the upper halo (see Section [6.1.2|) . then 
the interactions of the HVCs with the ambient gas may 
result in a large number of high ions at large distances 
from the plane. 

Our example composite model includes contributions 
from HVCs (this paper, specifica lly our referenc e model, 
Model B), extraplanar SNRs (iSheltonl 12006ft. radia- 
tively cooling Galactic fountains ( Shapiro fc Benjamin! 



tivery < 

11993!) . and photoion ization by an external radiation field 
(llto fc Ikeuchl Il988h . The contributions from these var- 
ious components to the observed column densities are 
presented in Table [3] (see the footnotes to that table for 
more details) . Note that the overall normalization of the 
Galactic fountain component (component 3 in Table [3]) 
is essentially a free parameter, as the mass flow rate in 
the fountain is not well known. We the r efore rescaled 
the predictions from IShapiro fc Benjamin! (|1993l ) so that 
this component would account for all of the O VI not ac- 
counted for by HVCs and SNRs (this rescaling amounted 
to multiplying their predictions by 0.24-0.38, depending 
on the particular model in their Table 1). Hence, our 
model accounts for all of the observed O VI by design. 

In our composite model, HVCs, SNRs, and Galac- 
tic fountains (components 1-3) can account for approxi- 
mately half of the observed C IV and N V. The Galactic 
fountain model of iShapiro fc Benjamin! (|1993f ) includes 
the effects of photoionization by the cooling gas's own 
radiation, but not by an external radiation field. Pho- 
toionization by the extragalactic background, or by ra- 
diation from hot stars escaping from the disk, may be 
able to account for the shortfall i n C I V. For example, 
the hybrid model of llto fc Ikeuchil (|198 8l includes a pho- 
toionized component with T = 10 4 K, which could ac- 
count for «40% of the observed C IV (component 4 in 
Table [3]); i.e., most or nearly all of the observed C IV is 
now accounted for. However, this photoionized compo- 
nent is unlikely to be able to account for the additional 
N V, due to N V's higher ionization potential. 

The shortfall in N V in our composite model is 
due to the fact that the various components of our 
composite model all underpredict the N V column 
density relative to the O VI column density. The 
predicted ion ratios for HV Cs (this paper), extra- 
planar SNRs (jSheltonl 12006( 1 . and Galactic fountains 
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Table 3 

Composite Model of Low- Velocity High Ions in Halo 





CIV 


N V 


o vi 


Observed 3 - (10 13 cm~ 2 ) 


10 


2.5 


15.5 


(1) HVCs b 


13% 


22% 


37% 


(2) Extraplanar SNRs c 


6% 


12% 


17% 


(3) Galactic fountains' 1 


23-37% 


14-21% 


46% 


(4) Photoionization by 








external radiation field 


39% 






Total 


81-95% 


47-55% 


100% f 



Note. — Each model component's contribution is ex- 
pressed as a percentage of the observed column density, in 
the first row. Any apparent discrepancies between the in- 
dividual components' contributions and the totals are due 
to rounding. 

a C IV: ISavage He Wakkerl d2009l') ; N V: ISavaee et~aLl 119971) ; 
O VI: ISavage et al.l l|2003l), after removin g contribution 
from the Local Bubble {Ocgcrlc et al. 200TJ. 
b Predictions from Model B (our reference model); specif- 
ically, the average of the Domain-only and Domain + Es- 
c aped prediction s from the lower half of Table [2] 
c ISheltonl {2005), Table 8. We have chosen the model with 
the median O VI prediction (case 1, drag coefficient 1). 
d Shapiro & Benjamin (1993), Table 1, rescaled such that 
this component accounts for all of the O VI not already 
accounted for by HVCs and SNRs. This model includes 
t he effects of se lf-photoionization. 

e llto He Ikeuchil ffegiD . Figure 3(a). We assume that all of 
the C IV in this figure is due to the photoionized component 
of their model, and that this component produces negligible 
N V and O VI. 

f Our composite model reproduces 100% of the O VI by 
design. 

(jEdgar & Chevalierl [19861: iShapiro fc Benjamin! [19931) 
are all in the range log [N(N V)/JV(0 VI)] « -1.3 
to —1.0. In comparison, the observed value is 
log[iV(N V)/iV(0 VI)] = -0.8 (see Tabled. Some 
additional source ol gas with T w 2 x 10 5 K would 
be needed to increase the N V column density rela- 
tive to the C IV and O VI values, although what this 
source could be is not obvious. Plane-parallel turbu- 
lent mixing layer models tend to give higher N V-to- 
O VI ratios (log [ V(N V)/iV(Q VI)] = - 0.9 to -0.4; 
iSlavin et all 119931: iKwak fc Sheltonl I2010D . but as the 
composite model already includes mixing between HVC 
and ambient gas, it is not clear how more turbulent mix- 
ing layers could be added to the model. 

7. SUMMARY 

We have presented further results from our NEI hydro- 
dynamical simulations of HVCs (Paper I). In this paper, 
we concentrated on the low-velocity high ions that re- 
sult when an HVC passes through the hot halo. These 
high ions arise from the mixing of cool cloud gas with 
hot ambient gas. Initially, this mixed high-ion-rich gas 
is located near the HVC, and travels with HVC-like ve- 
locities, but later it falls behind the HVC and slows to 
ISM-like velocities, while retaining its high ion content. 

We examined a suite of seven models, covering differ- 
ent cloud velocities, cloud densities, cloud density pro- 
files, and cloud sizes. In general, the cloud velocity and 
density profile have little effect on the masses of the 
high ions that result. Larger or denser clouds result in 
greater masses of high ions, as one would qualitatively ex- 
pect. An important result is that, except for the fastest 
model HVC (|u| = 300 km s -1 ), the quantities of low- 



velocity ions are generally larger than the quantities of 
high- velocity ions. This result suggests that HVCs could 
be an important source of low-velocity high ions in the 
halo. 

We examined this suggestion more quantitatively, us- 
ing the HVC infall rate to estimate the average column 
densities of low-velocity high ions in the halo due to 
HVCs. After accounting for the contribution from ions 
at times beyond the ends of our simulations, and being 
conservative regarding the HVC infall rate, we find that 
our models can account for >30% of the O VI column 
density observed in the halo. This implies that the colli- 
sionally ionized gas in material shed by HVCs is a signifi- 
cant source of the low- velocity O VI observed in the halo. 
In contrast, such gas is probably not a significant source 
of low-velocity C IV: our reference model can account 
for only 12-14% of the observed column density. This is 
probably because the observed C IV is likely affected by 
photoionization, which our model does not include. 

We used the predictions of our HVC model in a simple 
composite model of the low- velocity high ions in the halo, 
in which we combined the contributions from HVCs, ex- 
traplanar SNRs, radiatively cooling fountain gas, and 
photoionization from an external radiation field. By de- 
sign, this model accounted for all of the observed O VI. 
We found that the model could account for most or all of 
the observed C IV, but only about half of the observed 
N V. It is not obvious what the source of the additional 
N V could be. Although this composite model was con- 
structed in a relatively simple way, and failed to account 
for all the observed ions, we emphasize the point that 
any complete model of the high ions in the halo should 
include the contributions from HVCs, particularly to the 
column density of low-velocity O VI. 

The software used in this work was in part devel- 
oped by the DOE-supported ASC/ Alliance Center for 
Astrophysical Thermonuclear Flashes at the University 
of Chicago. The simulations were performed at the Re- 
search Computing Center (RCC) of the University of 
Georgia. We acknowledge use of t he R software package 
(|R Development Core Teamll2008l ). This work was sup- 
ported by NASA grant NNX09AD13G, awarded through 
the Astrophysics Theory and Fundamental Physics Pro- 
gram. DBH acknowledges funding from NASA grant 
NNX08AJ47G, awarded through the Astrophysics Data 
Analysis Program. 

APPENDIX 

A SIMPLE ANALYTICAL MODEL OF A SPHERICAL HVC 

Here we describe a simple analytical model of a spheri- 
cal HVC that loses mass at a rate proportional to its sur- 
face area (introduced in Section 3.3.3 of Paper I). This 
model is helpful for comparing the results from hydro- 
dynamical simulations of different-sized HVCs (see Sec- 
tions [321 and [4~2]) • 

We assume that the mass of low- velocity high ions of 
a given type, M- lon , increases with time at a rate that 
is proportional to the rate at which the HVC sheds 
its mass, and thus to the HVCs surface area, i.e., 
dMi on /dt = ar 2 (t). In this simple model, the cloud ra- 
dius decreases linearly with time from its initial value, ro 



16 



(see Paper I), and so we can rewrite the mass loss rate as 
dM ion /dt = a(r — kt) 2 . Note that the constants a and k 
are unimportant, as long as they are independent of r*o. 
Note also that we are ignoring the fact that the ions of 
interest will subsequently recombine or ionize (i.e., there 
should be an additional, negative term on the right-hand 
side of the expression for dMi on /dt). 

We integrate dM- lon /dt with respect to t, with the 
boundary condition M[ on = at t = 0. In order to com- 
pare different-sized clouds at equivalent stages of their 
evolution, we express the solution as a function of the 
rescaled time, t/ro, obtaining 



M ion (i/r ) = ^ 



1 



1 



(Al) 



Hence, we expect that 

M ion (t/r ) oc oc M H vc,o- 



(A2) 

We use this result in Section [3T2l Note that this result is 
only valid if t < ro /k and if the rate at which the ions of 
interest subsequently recombine or ionize is negligible. 

We can also use this model to compare the ion column 
densities expected from different-sized clouds. If the end 
time, tf, for the integral in Equation ([3]) is chosen such 
that tf/rg is the same for different-sized clouds, then we 
expect that 



TV(ion) oc 



Jo M ion dt 

Jvi HVC.O 



Jo tf/r ° M ion (t/r Q )d(t/r ) 



■■ r 



/u HI 

Jvi HVC.O 



ocr , 



(A3) 

where, for the final step, we have used the fact that 
Mi on (t/ro) oc Mjjy C0 in this simple model (Equa- 
tion (|A2[) ) . We use this result in Section T4.2I Again, 
note that this result is only valid if the rate at which the 
ions of interest ionize or recombine is negligible. 
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